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SUl-iMARY 


The  subsonic  ^eroelastlo  stability  of  a  two^ 
dimensional  panel  resting  on  a  continuous  elastic  foundation 
was  investigated  both  experimentally  and  theoretically. 

Experimentally,  tests  were  conducted  on  a  104  in. 

X  24  in.  X  0.018  in.  rectangular  aluminum  panel  in  a  low- 
speed  wind  tunnel.  Definite  flutter  of  a  travelling  wave 
type  was  observed.  Films  and  oscillograph  records  were 
taken. 

Theoretically,  a  finite  panel,  two  mode,  standing- 
wave  analysis  was  shown  to  give  essentially  the  same  be¬ 
havior  as  the  infinite  panel,  travelling-wave  analysis  of 
Miles  (Ref.l)  for  this  panel  on  an  elastic  foundation.  Al¬ 
though  a  mild,  divergence-type  instability  exists  for  these 
panels,  the  dominant  instability  was  shown  to  be  of  a  travel 
ling  wave,  flutter  type.  The  role  of  additional  internal 
damping  was  Investigated  and  found  to  have  an  interesting 
destabilizing  effect  between  the  divergence  and  flutter 
speeds. 

Comparison  of  experiment  and  theory  showed  good 
agreement  in  the  prediction  of  the  flutter  speed  and  wave¬ 
length,  but  rather  poor  agreement  in  the  wave  speed  and  fre¬ 
quency  at  flutter.  This  discrepancy  was  attributed  to  vari¬ 
ous  limitations  in  the  test  set-up  as  well  as  to  the  general 
difficulty  of  predicting  the  wave  speed  and  frequency  as  ac¬ 
curately  as  the  flutter  speed. 

The  present  investigation  should  be  of  interest  in 
problems  of  hydroelasticity,  axially  symmetric  cylinders, 
and  inflatable  structures  at  low  speeds,  as  well  as  to  the 
general  problem  of  panels  lying  on  springy  elastic  materials 
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1.  INTRODUCTION 

"  ^THe  problOT  considered  here  is  the  aeroelastlc 
stability  of  a  long  panel  resting  on  a  continuous  elastic 
foundation.  Such  a  situation  might  arise  for  a  panel, 
lying  on  some  springy  insulating  material,  exposed  to  an 
air  flow  over  its  top  side. 

The  theory  of  Miles  in  Ref.l  for  infinitely  long 
panels  can  readily  be  adapted  to  describe  this  behavior.  In 
this  theory,  unstable  travelling  waves  ara  predicted  above 
a  critical  flutter  speed  for  both  the  subsonic  and  super- 
sonic  cases.  For  panels  resting  on  an  elastic  foundation, 
these  unstable  travelling  waves  possess  a  finite  wavelength, 
and  hence  give  the  possibility  of  actually  viewing  these 
waves  on  a  long  panel.  Accordingly,  it  was  decided  to  build 
such  a  panel  and  test  it  In  a  conveniently  available  low- 
speed  wind  tunnel.  It  was  hoped  to  shed  some  light  on  the 
travelling  wave  versus  standing  wave  theories  of  panel  flut¬ 
ter,  and  also  to  reinvestigate  the  possibility  of  panel 
flutter  at  subsonic  speeds. 

The  present  report  gives  an  account  of  the  experi¬ 
ments  performed,  and  a  detailed  theoretical  analysis  of  the 
panel  from  both  the  infinite  length  and  finite  length  point 
of  view.  Some  earlier  preliminary  work  in  this  general 
area  was  done  by  Dowell  in  Ref. 2. 

The  results  of  the  present  Investigation  may  also 
be  of  interest  in  the  related  problems  of  hydroelastic  panel 
flutter,  axially  symmetric  cylinder  flutter  at  low  speeds, 
and  panel  flutter  of  inflatable  structures  at  low  speeds. 
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2.  EXPERIMENT 


2»1  Panel  Test  Model 

The  model  Co  be  tested  consisted  basically  of  a 
flat  rectangular  sheet  of  aluminum  2024-T3  with  dimensions 
104  in.  X  24  in.  x  0.018  in.  The  panel  was  orthogonally 
stiffened  by  a  series  of  aluminum  angle  stiffeners,  spaced 
1.58  inches  apart,  and  running  perpendicular  to  Che  long 
edges  of  Che  panel.  The  panel  was  attached  to  a  rigid 
wooden  box  frame  by  resting  on  three  rows  of  helical  springs 
running  parallel  Co  Che  long  edges  of  the  panel,  and  by  a 
plastic  hinge  along  one  of  the  short  sides  (the  leading  edge 
of  Che  panel) .  The  two  long  edges  of  Che  panel  were  left 
free,  while  the  other  short  side  (the  trailing  edge)  was 
left  free  during  some  of  Che  tests,  and  pinned  during  others. 
A  sketch  of  Che  panel  and  its  support  system  is  shown  in 
Ftg.l. 

Both  Che  aluminum  angle  stiffeners  and  Che  helical 
coil  springs  were  attached  to  Che  aluminum  sheet  by  Armstrong 
A-2  epoxy  resin  glue.  This  provided  a  strong  and  rugged  bond 
The  angle  stiffeners  were  used  in  order  to  insure  txra-dimcn- 
sional  behavior  of  the  panel.  The  springs  in  each  of  Che 
three  rows  were  spaced  1.58  Inches  apart  in  the  stream  di¬ 
rection.  Since  the  wavelengths  anticipated  were  of  Che 

During  the  early  part  of  the  test  program  (Run  17),  the 
very  rear  part  of  the  panel  was  damaged.  The  panel  was  then 
shortened  from  104  in.  to  93  in.  for  the  remainder  of  the 
tests. 

An  oil-can  effect  existed  originally  in  the  forward  part 
of  Che  sheet.  This  necessitated  somewhat  larger  size  stif¬ 
feners  in  this  vicinity  in  order  Co  get  rid  of  this  three- 
dimensional  effect. 
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order  of  24  Inches^  it  was  felt  that  the  discrete  springs 
adequately  represented  a  continuous  elastic  foundation. 

The  leading  edge  was  pinned  by  a  mylar  strip 
hinge.  The  trailing  edge,  which  was  kept  free  for  much  of 
the  test  program,  was  limited  to  about  1/2  in.  vertical 
motion  above  the  undeformed  position  by  three  loosely  fit¬ 
ting  bolts.  These  were  placed  In  after  Run  17,  when  ex¬ 
cessive  motion  damaged  the  trailing  edge.  During  the  latter 
part  of  the  test  program,  the  trailing  edge  was  pinned  In¬ 
stead  of  free.  This  was  provided  by  replacing  the  three 
loosely  fitting  bolts  by  three  plastic  knife  edges.  Slotted 
holes  were  made  In  the  panel  to  prevent  axial  tension  forces 
from  developing  when  the  panel  deflected  vertically. 

The  bottom  of  the  wooden  box  frame  was  covered  by 
a  wooden  panel  with  numerous  one-inch  diameter  holes  to  al¬ 
low  tvmnel  static  pressure  to  be  developed  on  the  Inside  of 
the  aluminum  test  panel. 

A  square  grid  of  black  lines  was  painted  on  the 
panel,  for  ease  in  viewing  and  observing  the  panel  behavior. 
The  distance  between  lines  was  roughly  six  inches,  and  the 
lines  perpendicular  to  the  flow  were  numbered  from  0  to  16, 
starting  at  the  leading  edge.  Eight  SR-4  strain  gages  were 
also  placed  on  the  underside  of  the  panel  near  its  mid- 
section.  These  were  to  be  used  to  obtain  frequency  and 
phasing  behavior.  The  location  of  the  grid  and  the  strain 
gages  is  shown  in  Fig. 2. 

2.2  Static  and  Vibration  Tests 

Prior  to  flutter  testing  in  the  tunnel,  static 
tests  were  made  on  a  typical  helical  coil  spring,  and  vibra¬ 
tion  tests  were  made  of  the  panel  as  a  whole  in  its  wooden 
box  frame. 
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The  helical  coil  springs  used  had  a  length  of  11/2 
inches^  a  dianeter  of  1/2  inch,  and  could  be  compressed  to 
bottom  without  buckling.  The  results  of  a  static  compres¬ 
sion  and  tension  test  of  a  typical  spring  is  shown  in  Fig. 3. 
The  spring  is  seen  to  be  fairly  linear  over  a  wide  range 
and  has  a  spring  rate  of  0.88  Ib/in. 

The  panel  was  vibrated  in  its  wooden  supporting 
frame,  by  attaching  an  elcctrodynamic  shaker  to  the  panel 
centerline,  about  8  inches  from  the  leading  edge.  The  natur¬ 
al  modes  and  frequencies  were  obtained  by  visual  observation. 
Salt  was  sprinkled  over  the  panel  to  determine  the  position 
of  the  node  lines,  and  the  audio  oscillator  output  frequency 
was  checked  with  a  strobotac.  The  experimental  modes  and 
frequencies  together  with  theoretical  ones  computed  from  the 
spring  and  mass  properties  (see  Section  3.1),  are  given  In 
Fig. A  for  the  higher  modes  (n  ■»  7  to  n  “  11).  The  observed 
frequencies  are  in  reasonable  agreement  with  those  calculated. 
The  observed  node  lines  are  not  as  clear  cut  or  uniformly 
spaced  as  the  theoretical  ones.  However,  they  do  resemble 
the  calculated  ones,  and  the  trends  with  n  arc  quite  good. 

As  can  be  expected  from  the  theoretical  analysis,  the  lower 
modes  (n  ■  1,  2,  3,  etc.)  are  too  close  together  to  be  sepa¬ 
rately  excited  and  observed  by  this  single  shaker. 

The  panel  was  also  vibrated  with  the  shaker  located 
16.5  inches  from  the  trailing  edge.  As  these  mode  shapes  and 
frequencies  were  quite  similar  to  those  obtained  with  the 
shaker  in  the  forward  position,  they  are  not  repeated  here. 

2.3  Mounting  in  Tunnel 

The  test  panel  and  its  supporting  wooden  box  frame 
were  mounted  in  the  M.I.T.  Wright  Brothers  low-speed  wind 
tunnel.  This  tunnel  has  a  7  1/2  ft  x  10  ft  elliptical  test 
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section  which  Is  vented  to  approximately  atmospheric  pressure, 
and  is  capable  of  running  to  speed’s  of  about  140  mph  at  at¬ 
mospheric  pressure.  The  wooden  box^  frame  on  which  the  test 
panel  lay,  was  placed  between  two  end  plates  extending  from 
the  ceiling  to  about  three  feet  from  the  floor,  as  shovm  in 
Fig. 5.  These  end  plates  were  used  to  try  to  assure  two  di¬ 
mensionality  of  the  flow.  They  were  made  of  1/2-inch  ply¬ 
wood,  except  for  the  center  portion  of  one  side  which  was 
plexiglas  in  order  to  allow  visual  observation  of  the  test 
panel  surface  during  the  flutter  run.  The  wooden  box  frame 
could  be  changed  in  angle  of  attack  somewhat  in  order  to  ad¬ 
just  the  pressure  and  velocity  distribution  over  the  panel. 
Three  pitot-static  tubes  were  placed  about  12  inches  over 
the  centerline  of  the  panel  in  order  to  measure  the  velocity. 
Later  on  in  the  test  program,  nine  more  pitot-static  tubes 
were  added  in  order  to  survey  the  pressure  and  velocity  dis¬ 
tribution  over  the  panel  more  accurately. 

For  observing  visually  the  static  vertical  deflec¬ 
tions  of  the  panel  under  the  action  of  the  air  forces,  a 
series  of  six  long  black  lines  1/4-inch  apart  were  ruled  on 
the  inside  face  of  one  of  the  wood  beams  of  the  supporting 
wooden  box  frame.  The  static  vertical  deflection  of  the 
aluminum  panel  could  then  be  estimated  visually  by  reading 
it  against  this  grid. 

The  panel  surface  with  its  six-inch  square  grid 
could  be  viewed  visually  through  the  plexiglas  end  plate 
and  the  wind  tunnel  window.  Movies  were  taken  by  a  16  mm 
camera  from  this  position.  Also,  the  eight  strain  gages 
on  the  panel  underside  were  connected  to  a  recording  oscillo¬ 
graph. 
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2.4  Flutter  Test  Procedure 


During  a  typical  flutter  test  rxin,  the  tunnel 
speed  was  raised  to  a  given  value  and  held.  For  this 
speedy  the  tunnel  pitot:-’ static  tube  was  read.  This  gave 
the  reference  tunnel  setting  and  represented  the  dynamic 
pressure  in  inches  of  alcohol  at  some  reference  point  in 
the  tunnel  test  section.  The  three  pitot-static  tubes 
over  the  panel  itself  were  then  read^  in  order  to  obtain 
the  actual  velocity  at  the  panel.  The  static  deflection 
shape  of  the  panel  was  then  recorded  by  visual  observation 
of  the  deflections  of  the  long  edge  of  the  panel  against 
the  ruled  1/4-inch  grid  on  the  inside  face  of  the  supporting 
wooden  frame.  The  above  test  runs  were  repeated  in  small 
tunnel  velocity  increments  until  significant  vibrational 
activity  or  flutter  was  observed  visually  over  the  panel 
surface.  Then^  both  oscillograph  records  of  the  strain 
gagcs^  and  16  mm  movies  of  the  panel  surface  were  taken. 

A  test  log  of  the  different  runs  carried  out  Is 
given  In  Fig. 6.  This  test  program  is  a  consequence  of  at- 
temps  to  repeat  results^  to  Investigate  both  free  and 
pinned  trailing-edge  conditions,  to  minimize  the  effects  of 
nonuniform  pressure  and  velocity  distributions  over  the 
panel,  and  to  determine  the  detailed  pressure  and  velocity 
distributions  over  the  panel. 

During  Run  17,  the  panel  trailing  edge  fluttered 
violently  and  as  there  were  no  stops  at  that  time,  the  last 
nine  Inches  of  the  panel  was  damaged.  This  was  subsequently 
repaired  by  cutting  off  the  last  eleven  Inches  of  the  panel 
and  Installing  loosely  fitting  bolt  stops  (as  mentioned  In 
Section  2a) .  For  the  remainder  of  the  tests,  the  flutter 
that  did  occur  remained  nondestructive. 
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Early  in  the  test  program,  a  significant  static 
depression  occurred  near  the  leading  edge  of  the  panel  (see 
Fig. 9).  This  appeared  to  be  attributed  to  nonuniform  pres¬ 
sure  and  velocity  distributions  over  the  panel,  and  was 
remedied  somewhat  by  increasing  the  angle  of  attack  of  the 
panel  between  the  end  plates.  Later  on  in  the  test  program 
(Run  100  and  on),  more  accurate  pressure  and  velocity  sur¬ 
veys  over  the  panel  were  performed  using  twelve  pitot- 
static  tubes  in  place  of  the  original  three. 

2.5  Test  Results 

a)  Pressure  and  Velocity  Surveys 

The  results  of  pressure  and  velocity  surveys 
taken  during  Runs  100  through  125  are  shown  in  Figs. 7  and  8. 
These  were  obtained  for  the  different  angle-of-attack  set¬ 
tings  and  tunnel  settings.  There  was  some  small  scatter  in 
the  original  data  which  was  traced  subsequently  to  slight 
leaks  in  the  pressure  lines  of  a  few  of  the  pitot-static 
tubes.  The  data  in  Figs. 7  and  8  represent  the  best  faired 
estimate  of  the  pressiire  and  velocity  distributions  taken 
six  inches  above  the  surface  of  the  panel.  There  was  no 
significant  difference  between  the  distributions  six  inches 
above  the  panel  and  those  twelve  inches  above  the  panel. 

The  undesirable  variation  of  the  pressure  and  velocity  dis¬ 
tributions  over  the  panel  surface  are  to  be  noted.  On  the 
underside  of  the  panel,  the  pressure  distribution  was  found 
to  be  fairly  uniform  over  the  entire  length  of  the  panel  for 
each  given  angle-of-attack  setting  and  tunnel  setting. 

b)  Static  Deflections 

The  static  deflections  of  the  panel  for  different 
angle-of-attack  settings  and  tunnel  settings  are  shown  in 
Fig. 9.  The  prominent  depression  near  the  leading  edge  is 
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clearly  visible.  This  tended  to  become  smaller  for  the 
high  angle-of-attack  setting.  Also  clearly  visible  is  a 
sinusoidal  wave  pattern  of  about  two-ft  wavelength  that  . 
developed  near  the  trailing  edge  on  many  of  the  runs,  par¬ 
ticularly  on  Runs  114  through  121,  the  low  angle  of  attack, 
pinned  case. 

c)  Flutter  Observations 

Flutter  was  observed  for  many  of  the  runs.  Out¬ 
side  of  Run  17,  it  tended  to  be  nondestructive  and  limited 
in  amplitude,  so  that  it  could  be  vif'wed  for  relatively 
long  periods  of  time.  The  double  amplitude  reached  values 
of  one-half  inch  on  many  of  these  runs.  The  occurrence  of 
flutter  is  indicated  in  the  test  log.  Fig. 6.  These  results 
are  summarized  lie  Fig. 10.  The  flutter  velocity  quoted  was 
taken  at  70  inches  aft  of  the  leading  edge  (roughly  at  3/4 
of  the  length  of  the  panel) . 

Considerable  variation  In  the  flutter  speeds  oc¬ 
curred  particularly  between  the  different  angle-of-attack 
settings.  There  was  a  lesser  change  between  the  free  and 
pinned-end  conditions.  It  is  felt  that  the  results  for  the 
high  angle-of-attack  case  are  the  most  reliable  since  these 
had  the  smallest  static  depression  at  the  leading  edge,  and 
the  smallest  pressure  and  velocity  variation  over  the  panel 
surface.  The  large  static  depression  in  addition  to  aero¬ 
dynamic  effects,  shortened  the  effective  length  of  the  panel 
and  may  have  introduced  midplane  tension  and  structural  non¬ 
linear  effects  into  the  panel  behavior. 

d)  Movies 

Movies  of  many  of  the  flutter  test  runs  were  taken 
with  a  16  mm  camera  at  64  fraroes/sec.  The  more  interesting 
shots  were  compiled  together  into  a  single  400-ft  film  strip 
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representing  22  different  runs.  The  film  log  for  these  pic¬ 
tures  is  given  in  Fig. 11.  The  movie  films  recorded  both 
steady  flutter,  and  the  onset  of  flutter  in  some  cases. 
Typical  flutter  film  sequences  are  shown  in  Figs. 12a  through 
12d  for  the  high  and  medium  angle-of-attack  cases.  An  ex¬ 
ample  of  the  static,  divergence-type  phenomenon  observed 
during  Runs  114  through  121  is  given  in  Fig.l2e. 

Analysis  of  these  flutter  film  sequences  and  other 
runs  in  the  film  strip  seem -to  indicate  a  flutter  frequency 
of  about  14  cps,  a  wavelength  of  about  2  to  2  1/2  feet,  and 
a  wave  speed  of  about  35  ft/sec.  The  waves  generally  ap¬ 
pear  to  be  travelling  downstream,  although  there  seems  to 
be  some  small  "sticking"  or  standing  wave  component  present 
in  some  of  them.  Some  small  distortion  of  the  travelling 
waves  from  a  pure  sinusoidal  shape  is  evident,  particularly 
for  the  larger  anqplitudes  of  wave  motion.  A  camera  speed 
somewhat  faster  than  the  64  frames/sec  used  here,  would  have 
helped  define  the  wave  motion  more  accurately. 

e)  Oscillograph  Records 

Oscillograph  records  of  many  of  the  flutter  test 
runs  were  taken.  The  oscillograph  log  is  given  together 
with  the  film  log  in  Fig. 11.  Typical  oscillograph  records 
for  the  same  runs  as  the  film  strips  of  Figs. 12a  Co  12d  arc 
given  in  Figs. 13a  through  13d.  On  these,  the  eight  strain- 
gage  traces  are  shown,  giving  clearly  the  frequency  and 
phasing  of  Che  different  locations.  Since  Che  gages  are  all 
roughly  similar  in  sensitivity,  a  measure  of  the  amplitudes 
is  also  obtained  for  these  locations.  To  compare  amplitudes 
between  different  oscillograph  records,  all  amplitudes  on  a 
given  record  should  be  multiplied  by  the  amplitude  factor 
A.F.  indicated,  in  order  to  adjust  roughly  for  the  different 
amplifier  gain  settings  used  during  Che  test  program. 


The  traces  shown  In  Figs* 13a  through  13d  show 
strong  evidences  of  standing  waves  with  some  small  tr^el- 
ling-v/ave  components  present.  The  frequencies  average  14,4 
and  13.7  cps  for  the- high  and  medium  angle-of -attack  settings, 
respectively.  The  half  wavelengths  average  about  four-gage 
lengths  or  12.8  inches.  This  is  particularly  clear  from 
Fig. 13d. 

Two  other  oscillograph  records,  from  Runs  12  and 
76,  are  shown  in  Figs.l3e  and  13f.  These  were  taken  slightly 
below  the  visually  observed  flutter,  and  show  here  clear  evi¬ 
dences  of  travelling  waves. 

f )  Summary  of  Results 

It  appeared  that  definite  flutter  was  obtained  for 
this  panel  at  speeds  of  about  95  mph,  frequencies  of  14.4  cps, 
wavelengths  of  about  2  to  2  1/2  feet,  and  wave  speeds  of  about 
35  ft/sec.  The  waves  generally  appeared  to  be  travelling 
downstream  although  there  seemed  to  be  some  "sticking"  or 
standing  wave  components  present  near  the  center  of  the  panel. 
The  flutter  speed  was  not  affected  much  by  the  tralllng-edge 
condition,  i.e.,  whether  pinned  or  free  there.  Prior  to 
flutter,  the  panel  assumed  a  sinusoidal  static  deflection 
shape  of  about  two  feet  wavelength  in  the  vicinity  of  the 
trailing  edge  (in  one  case.  Run  121,  this  became  a  rather 
sizable  deflection) .  There  also  appeared  a  prominent  static 
depression  in  the  vicinity  of  the  leading  edge  which  was 
remedied  somewhat  by  altering  the  angle  of  attack  and  hence 
pressure  distribution  over  the  panel. 

2.6  Improvement  of  Experimental  Test 

After  having  carried  out  the  above  experimental 
test  program,  it  has  become  apparent  that  several  improve¬ 
ments  can  be  made  to  obtain  more  accurate  test  data  in  any 
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future  work  along  these  lines. 

First,  some^/hat  larger  angle  stiffeners  should  be 
used  to  insure  two-dimensional  behavior  of  the  panel  and  to 
eliminate  any  possible  oil-canning  effects  in  the  panel  be¬ 
havior. 

Second,  more  care  should  be  given  to  maintaining 
uniform  pressure  and  velocity  distributions  over  the  panel, 
and  to  eliminating  the  static  depression  near  the  leading 
edge.  This  could  be  accomplished  by  more  anglc-of-attack 
manipulation,  a  longer  nosepiece  on  the  wooden  box  frame, 
bettor  contouring  of  the  walls  above  and  alongside  the 
panel,  and  better  venting  between  the  upper  and  lower  sides 
of  the  panel.  A  thorough  pressure  and  velocity  survey 
should  be  made  before  passing  on  to  the  flutter  tests. 

Third,  better  observation  techniques  could  be 
utilized.  These  would  include  using  a  faster  frame  speed 
camera,  better  calibration  of  the  strain  gages,  accurate 
synchronization  of  camera  with  oscillograph  records,  and 
photographing  more  of  the  panel  surface  in  the  movies. 
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3.  THEOIY 

3.1  .Free  Vibration  Modes 

The  equation  for  the  derlection  w(x,t)  of  a  two- 
dincnsional  panel  resting  on  a  continuous  clastic  founda¬ 
tion  iSj 

^  Kw  ==  Ap-  -  na 


where  for  the  present  panels 


D  =•  “  flexural  rigidity  =*  5.61  lb-in 

K  =»  foundation  stiffness  =  0.0696  Ib/in® 

tn  =  mass  per  unit  area  =*  5.79  x  10  *  lb-sec*/in® 

(2) 

£.p  »  pressure  loading  =  Ibs/in® 

I  “  length  of  panel  =*  104  in 

»  air  density  at  sea  level  •  0.1148  x  10  ®  Ib-scc^/ir"* 

For  free  vibrations  in  a  vacuum,  the  pressure  loading  is  as¬ 
sumed  zero.  Upon  utilizing  the  standard  separation  of  the 
variables  technique,  one  arrives  at  the  general  solution, 

W(x,i)  —  A  (cot  ♦  i)  f  C  *  3)  E  cara  M.X. 

■  (3) 


where 


U.  = 


mis’- K 


i*utting  in  the  boundary  conditions. 


@X  =  0  j'wrO 


one  obtains  the  characteristic  equation, 

LciL  ”  ”tiVr>-2^ 


Solving  the  above  for  ul,  and  making  use  of  Eq.U, 
following  natural  frequencies  are  obtained. 


0) 


mjt* 


the 


(5) 


The  corresponding  natural  mode  shapes  are. 


a) 


For  the  higher  values  of  n,  these  modes  resemble  pure  sine 
waves  except  in  the  immediate  vicinity  of  the  trailing  edge. 
The  n^b  mode  has  n  node  lines  spaced  approximately  a  distance 
J^n  +  1/A)  apart,  starting  from  the  leading  edge. 

For  the  present  panel,  the  theoretical  frequencies 
and  node  lines  for  free  vibrations  in  a  vacuum  are  shown  in 
Fig. A,  corresponding  to  values  of  n  ■  7  to  n  *»  11.  The  fun¬ 
damental  frequency  for  n  »  1  is  calculated  at  17.3  cps.  It 
is  seen  that  the  lower  modes  are  spaced  very  close  together 
in  frequency. 


For  the  actual  vibrations  of  the  panel  in  still 
air,  an  additional  virtual  mass  of  air  should  be  included 
with  the  panel  mass.  This  is  shown  in  Section  3.3  to  re- 
suit  in  a  total  mass  of  approximately. 


These  still  air  frequencies  are  also  shoxrm  in  Fig. A.  The 
mode  shapes  are  unchanged  by  this  additional  air  mas. 
. — 

For  still  air  on  both  sides  of  the  panel.  Also  for  high 
n,  X'»horo  sin  mrx/l  is  a  reasonable  approximation  to  the 
mode  shape. 
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3.2  .  I.nf inij:_e_ ^.npl  In  Incompressible  Flow 
a)  No  Internal  Damping 

The  flutter  of  an  infinite  panel  in  two-dimensional 
incompressible  flow  has  been  investigated  by  Milos  in  Ref.l. 
Ho  utilizes  there  an  axis  system  which  is  fixed  to  the  air 
at  rcst^  and  considers  the  infinite  panel  to  be  moving  with 
velocity  U  in  the  negative  x  direction.  The  equation  of 
motion  for  such  an  undamped  panel  on  an  clastic  foundation 
is, 

-  Kw  ,  =  Ap. 

Miles  has  shown  that  for  travelling  wave  solutions 
of  the  form, 

=  w,  e 

the  aerodynamic  pressure  distribution  for  air  acting  on  the 
top  side  of  the  panel  only,  is 

^  Iffl(ct-X) 

Ap  =  f»  ^  K  a  (11) 

Substituting  Eqs.lO  and  11  into  Eq.9,  one  is  lead  to  the 
equation 

f?”  l)"  -  =  °  (12) 

where 

The  pressure  acting  on  the  bottom  side  of  the  panel  is  as¬ 
sumed  constant  at  p^  in  the  present  analysis.  The  presence 
of  air  on  the  bottom  side  of  the  panel  is  examined  in 
Appendix  A  and  is  found  to  be  small  for  the  low  M.'s  of  in¬ 
terest  hero. 


lA 


?  =  F 


Co 

u 

Mr  m 


c_  = 


-  f  mtV 

m  V  ^  y 


(13) 

(14) 

(15) 

(16) 


In  the  above  quantities,  c  represents  the  wave  speed  with 
respect  to  the  air  at  rest,  c©  represents  the  wave  speed 
in  the  absence  of  air  (U  *•  0,  P*,  *0),  and  u  represents 
the  nass  density  ratio  of  the  panel  (for  air  on  one  side 
of  the  panel  only).  Equation  12  is  solved  for  Any 
solutions  for  ^  having  negative  imaginary  parts  represent 
instabilities. 


Solving  Eq.l2  explicitly  for  C>  Miles  obtains  the 
following  expressions  for  the  wave  speed  relative  to  still 
air,  c. 


1 


(17) 


and  for  the  wave  speed  relative  to  the  panel,  c  +  U, 

c.v  =  j  (18) 

It  is  readily  seen  from  these  expressions  that  negative  Imaginary 
values  of  c,  hence  instability,  will  occur  for 


Note  this  definition  of  ^  is  the  reciprocal  of  another 
definition  of  mass  density  ratio  frequently  used  in  flutter 
work. 
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The  above  criterion  is  plotted  in  Fig. 14,  The  corresponding 
v;avo  speed  relative  to  the  panel  is. 


c  >  U  = 


I  +  Y" 


In  applying  the  above  criterion,  it  is  Co  be 
noted  that  both  p.  and  Cq  arc  functions  of  the  wavelength  X. 
Hence,  one  must  investigate  the  variation  of  flutter  speed 
with  wavelength.  For  a  panel  on  an  elastic  foundation,  the 
variation  of  Cq  with  X  is  given  by  Eq.l6  and  is  shown  plotted 
in  Fig. 15.  The  corresponding  variation  of  flutter  speed  Up 
with  X  is  shown  in  Fig.  16.  Unlike  the  case  of  a  membrane  or 
a  plate  where  Che  minimum  flutter  speed  occurs  for  X-*  ©o, 
the  panel  on  an  clastic  foundation  possesses  a  minimum 
flutter  speed  at  a  finite  value  of  X  and  thereby  gives  hopes 
of  actually  seeing  these  travelling  waves  on  a  long  panel. 
This  wavelength  X  for  the  minimum  flutter  speed  is  given  by 

f,  (21) 

where  fi  is  a  correction  factor  given  by 


- 


and  is  seen  to  vary  from  fj  »•  1.0  at  high  values  of  p  to 
fi  ■  1.316  as  u  -►  0. 

For  the  present  panel,  the  characteristics  at 
flutter  are  calculated  numerically  as 
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Flutter  speed  Up  *  105.5  mph 

Wavelength  X  ■  25  in 

Wave  speed  c  +  U  -  11.4  ft/sec  (23) 

A 

Frequency  to  “  (c  +  U)  «  34.4  rad/sec  -  5.5  cps 

Density  ratio  ji  *  0.079 

Below  the  flutter  speedy  a  multitude  of  undamped  travelling 
wave  solutions  are  possible  corresponding  to  different  wave¬ 
lengths  X.  The  wave  speeds  relative  to  the  panel,  c  +  U, 
are  plotted  in  Fig. 17  for  different  wavelengths.  These 
were  computed  from  Eq.l8. 

It  is  of  Interest  to  investigate  the  unstable  be¬ 
havior  of  the  panel  above  the  flutter  speed  and  to  determine 
how  fast  the  instability  sets  in.  In  terms  of  a  coordinate 
system  X  fixed  to  the  moving  panel,  the  travelling  wave  solu¬ 
tions  of  Eq.lO  can  be  written 

wfM)  =  KC  e 

where  Cj^  and  Cj  are  the  real  and  imaginary  parts  respectively 
of  c.  The  quantity  (cj^  +  U)  ■  (c  +  U)j^  represents  the  wave 
speed  relative  to  the  panel,  while  the  quantity  Cj  ■  (c  +  U)j 
is  related  to  the  rate  of  growth  of  the  waves  as  they  move 
along  the  panel.  As  a  quantitative  measure  of  this  rate  of 
growth,  the  amplification  A  is  introduced.  This  amplifica¬ 
tion  A  is  defined  as  the  ratio  of  the  wave  amplitudes  during 
the  time  it  takes  the  wave  to  travel  one  wavelength  X  (i.c., 

A  2 wave  doubles,  etc.).  From  Eq.24  it  can  be  deduced 
that. 
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For  this  panel,  the  complex  wave  speed  c  +  U  is 
shoxm  plotted  in  the  complex  plane  in  Fig. 18  for  various 
speeds  and  wavelengths  X.  The  corresponding  amplification 
A  is  given  in  Fig. 19  also  for  various  speeds  and  wavelengths. 
It  can  be  seen  that  the  instability  sets  in  very  sharply 
after  Up  is  exceeded,  and  the  amplification  A  reaches  values 
of  the  order  of  100  for  2  or  3  mph  above  the  flutter  speed. 


The  preceding  analysis  was  for  a  panel  with  no  in¬ 
ternal  damping  present  in  the  structure.  Next,  consider  the 
effect  of  adding  some  internal  damping  to  the  structure. 

The  basic  equation  of  motion  would  then  be, 

and  the  corresponding  equation  in  C  would  read 

(C'O  -  C  +(*■?-  “  °  <28) 

Solving  the  above  equation  for  C  i^^w  gives 
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('♦r)?  = -I- 4 («) 

where  C  “  BX/27raU  is  a  measure  of  Che  additional  damping 
present  in  the  structure.  Unstable  solutions  will  now  oc¬ 
cur  if  the  imaginary  part  of  the  radical  has  an  absolute 
value  greater  than  G/2. 

The  imaginary  and  real  parts  of  the  above  radical 
can  be  evaluated  as 


I  III.  I  I  I  I  I 

(30) 


R^{r] 


2-  xU\r] 


(31) 


Unstable  solutions  will  appear  for 


£  <  J- 


By  routine  algebraic  manipulation,  the  above  inequality  re¬ 
duces  simply  to. 


<  r 


(33) 


2 

At  the  value  Co  "  it  can  be  shown  by  direct  substitution 
that  Im  "  +  G/2.  It  is  interesting  to  note  that  the 

above  criterion  Zq.33  does  not  depend  on  the  magnitude  of 
the  damping  G  present  in  the  structure.  Hence  for  any  small 
amount  of  damping  present  instability  will  set  in  for 
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rather  than  the  somewhat  higher  value  given  previously  by 
Zq.l9  for  the  completely  undamped  case.  This  new  criterion 
Is  shown  plotted  in  Fig. 14  together  with  the  previous  criterion 
for  the  undamped  panel.  The  difference  between  the  two  crite¬ 
ria  increases  for  increasing  p.  At  the  onset  of  this  new  in¬ 
stability,  the  ”e  [  I  and  one  obtains  from  Fq.29 

the  wave  speed  relative  to  the  panel  c  +  U  as, 

c  +  U  =«  0  (35) 

This  indicates  a  static  divergence- type  instability  rather 
than  a  travelling  wave  flutter-type  instability  as  obtained 
in  the  undamped  case. 

The  presence  then,  of  even  a  very  small  amount  of 
internal  damping  Is  sufficient  to  lower  the  Instability 
speed  of  the  panel  and  to  change  its  character  from  a  travel¬ 
ling  wave  flutter-type  instability  to  a  static  divergence- 
type  instability.  This  behavior  has  been  noted  also  by 
other  authors,  for  example,  Leonard  and  Hedgepeth  in  Ilef.3 
for  the  similar  case  of  thin  circular  cylinders  in  an  air- 
stream,  and  also  more  recently  by  Landahl  in  Ref .4  in  a 
study  of  the  boundary  layer  on  a  flexible  panel. 

To  apply  Che  criterion  given  by  Eq.34,  one  must  in¬ 
vestigate  the  variation  of  Che  Instability  speed  with  wave¬ 
length  X,  as  was  done  for  the  undamped  case.  This  diver¬ 
gence  speed  is  shown  plotted  in  Fig. 16  together  with  the 
previous  flutter  speed  for  the  undamped  panel.  The  v7avo- 
length  X  for  the  minimum  divergence  speed  is  given  by  Eq.21 
with  the  correction  factor  fj  ■  1.316. 


For  the  present  panel  with  a  slight  amount  of 
damping  present,  the  characteristics  at  instability  .(di¬ 
vergence  now)  are  calculated  numerically  as 

Divergence  speed  Up  “  101.5  mph 

Wavelength  A  *  25  in 

Wave  speed  c  -i-  U  =*  0  (36) 

Frequency  m  **  0 
Density  ratio  p-  “  0.079 

The  above  divergence  speed  is  somewhat  lower  than  the  previ¬ 
ous  undamped  panel  flutter  speed  of  105.5  mph.  Below  the 
above  divergence  speed,  the  travelling  wave  solutions  cor¬ 
responding  to  different  wavelengths  arc  all  damped  travel¬ 
ling  waves. 

It  is  again  of  interest  to  investigate  the  un¬ 
stable  behavior  of  Che  panel  above  Che  divergence  speed 
and  Co  deCermine  how  fast  this  insCabiliCy  sets  in.  The 
amount  of  internal  damping  present  in  the  panel  foundation 
can  be  more  conveniently  expressed  as 

B  =  2.  (t>.R.)|k  rrT  (37) 


where  (D.R.)  represents  the  familiar  critical  damping  ratio 
and  is  of  Che  order  of  0.01  for  this  metallic  structure. 
Using  Eq.37,  Che  damping  parameter  G  of  Eqs.29  Co  32  be¬ 


comes 

G  =  x  nr 

IT  U  >1 


(38) 


For  the  present  panel,  Che  complex  wave  speed 
c  U  was  computed  from  Eqs.29  to  32  for  speeds  above  dl- 
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vergcnce,  and  for  several  values  of  damping  ratio  (D.R. ). 
The  results  arc  shown  plotted  in  the  complex  plane  in  Fig. 

20.  The  corresponding  amplification  A  is  given  in  Fig. 21. 
Only  the  behavior  of  the  minimum  wavelength  X  ■»  25  inches 
is  shown  since  the  other  X‘s  show  similar,  but  less  criti¬ 
cal  trends. 

From  Figs. 20  and  21,  the  transition  between  the 
undamped  and  the  damped  panel  now  becomes  clear.  It  is 
seen  that  for  low  damping  ratios,  the  divergence-type  in¬ 
stability  represented  by  Eq.34  is  a  very  mild  one,  and  the 
amplification  A  does  not  become  significant  until  the  air 
speed  approaches  that  for  the  flutter-type  instability 
represented  by  £q.l9.  Correspondingly,  the  wave  speed 
(c  +  U)j^  changes  from  zero  to  the  value  for  the  flutter- 
type  instability  given  by  Eq.20.  For  large  damping  ratios, 
the  divergence- type  instability  does  become  significant. 

It  is  associated  with  very  low  speed  travelling  waves  and 
hence  would  resemble  divergence.  The  very  interesting 
destabilizing  effect  of  damping  in  the  speed  region  between 
101.5  mph  and  105.5  mph  is  to  be  noted.  Similar  behavior 
has  also  been  discussed  by  Landahl  in  Ref .4. 

To  conclude  this  section,  it  might  be  appropriate 
to  give  here  basic  formulas  for  the  occurrence  of  flutter 
and  divergence  for  these  infinite  panels  on  elastic  founda¬ 
tions.  Simple  algebraic  manipulation  of  Eqs.l5,  16,  19, 

21,  22,  and  34  gives  the  dynamic  pressure  at  flutter  as 


and  the  dynamic  pressure  at  divergence  as 
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f  = 


.ns  K""* 


In  the  expression  forlP  «•  Up,  the  appropriate  value  of  p.  Is 
found  through  consideration  of  Eqs.lS,  21,  and  22.  As  In 
many  other  types  of  flutter  Instability,  the  dynamic  pres¬ 
sure  at  flutter  becomes  Independent  of  mass  density  ratio 
for  heavy  panels  and  low  density  air  (p.-^  0). 

3.3  Finite  Panel  In  Incompressible  Flow 

It  Is  of  Interest  to  analyze  also  the  finite 
panel  and  compare  results  with  those  obtained  for  the  In¬ 
finite  panel. 

The  differential  equation  to  be  solved  Is  given 
by  Eq.l  subject  to  the  appropriate  boundary  conditions. 

For  simplicity,  the  panel  will  be  assumed  pinned  at  both 

it 

ends.  The  aerodynamic  pressure  for  air  acting  on  the  top 
aide  only,  can  be  expressed  as 

A|.  = 


ttJ  V  -at  ix;  ^ 


The  higher  modes  are  of  interest  in  this  analysis,  and  for 
these,  the  different  boundary  conditions  do  not  affect  the 
modes  and  frequencies  very  much  except  In  the  Immediate 
vicinity  of  the  trailing  edge. 

itit 

The  pressure  acting  on  the  bottom  side  of  the  panel  is  again 
assisned  constant  at  p  .  The  presence  of  still  air  on  the  bot 
tom  side  of  the  panel  Is  examined  In  Appendix  A  and  is  found  to 
be  small  for  the  low  pi’s  of  Interest  here. 
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where  w  is  the  deflection  of  the  panel.  (See,  for  example. 
Ref. 5).  As  the  general  solution  of  the  integro-differentlal 
equation  given  by  Eqs.l  and  41  is  extremely  difficult,  a' 
Galerkin-type  solution  will  be  used  by  assvmiing  modal  solu¬ 
tions  of  the  form, 

N 

W(5c,t)  =  ^  |,(t)  np  (42) 


These  satisfy  the  pinned  end  boundary  conditions  at  x  »»  0 
and  X  “  i.  Further,  the  aerod3mamic  pressure  will  be 
assumed  to  be  that  for  an  infinitely  long  wavy  wall.  This 
is  a  reasonable  assumption  except  in  the  immediate  vicinity 
of  the  leading  and  trailing  edges,  particularly  for  the 
higher  modes  (large  n)  of  interest  here.  For  such  an  in¬ 
finite  wavy  wall  of  shape  given  by  Eq.42,  the  aerodynamic 
pressure  Lp  of  Eq.41  reduces  simply  to  (sec  also  for 
example  Ref. 6), 


Placing  this  expression  together  with  Eq.42  into  the  basic 
panel  Eq.l,  and  applying  Galerkin's  method,  results  in  the 
scries  of  ordinary  differential  equations 
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(44) 


n,r  “  1,  2,  •  •  •  N 

The  sunnat: Ion  above  is  taken  over  all  the  r  terms  for  which 
n  +  r  is  an  odd  integer.  It  is  convenient  to  introduce  the 


following  nondimensional 

quantities. 

t  = 

>1 1"  I 

Q  = 

T> 

K 

KJi^ 

p 

(45) 

VUth  these  definitions,  the  equations  of  motion,  Eqs.44 
become  nondimens ionalized  to. 


[u£]  ^ 

^  it" 


[„V.  K  -  ^  0 

■  '■*  0  »  ^  \ 


Solutions  of  the  abov^  equations  arc  sought  in  the  standard 

s 

fashion  in  the  form  c  .In  physical  time,  the  solutions 
s  c 

are  given  by  e  where 
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Consider  now  the  couplJ.n£  between  the  n  and 

S 

n  +  1  nodes.  Upon  introducing  c  into  Eq.A6  one  obtains. 


[  S ]  [  5'{.  %£  ]  *  K  -(-OTTQ]]  I, 


=  0 


Setting  the  deteminant  of  the  above  equations  equal  to 
zero  yields  the  characteristic  equation  in  “s 

*  Aj.s''  ♦  Ao  =  O 

where 


A.  = 


A.= 


K  -  nTrQ^[(y>i-!)  ir^  K  “ 
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The  roots  ^  of  the  characteristic  equation,  Eq.A9, 
are  examined  as  Q  increases  from  zero.  Real  positive  roots 
represent  divergence  conditions,  while  complex  roots  with 
positive  real  parts  represent  flutter  instabilities. 

First  the  static  divergence  behavior  is  examined. 
This  occurs  for  Ag  **  0  v;hich  implies 


4-  K 
nTT 


(5i; 


For  the  present  panel,  K  “  1.46  x  10^,  and  the  quantity  Q 
versus  n  is  plotted  in  Fig. 22.  The  minimum  occurs  for 
n  “  8  and  the  corresponding  value  of  Q,  U,  and  wavelength 
at  divergence  are 

Qjj  -  73,570 

Ujj  “  101.5  mph  (52) 

X  *»  2j?/n  »  26  in 

These  values  agree  well  with  those  obtained  for  the  in¬ 
finite  panel  and  given  by  Eq.36. 

Next,  the  complete  dynamic  behavior  of  the  system 
is  examined.  Equation  49  can  be  solved  to  give  the  four 
roots  "s  as 

(53) 


Flutter  (dynamic  instability)  will  occur  only  if  the  inner 
Aa  —  4A4A0  becomes  an  imaginary  number,  that  is. 


if 


A.  -  4  A.  A.  <0  (54) 
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since  only  then  will  complex  roots  for  s  appear  with  posi¬ 
tive  real  parts.  Using  the  values  given  by  Eqs.50,  the 
above  criterion  for  flutter,  Eq.54  becoraes 


I  TT 


(55) 


Since  the  left-hand  side  of  the  above  criterion  is  always 
positive,  it  is  obvious  that  flutter  can  only  occur  for 
negative  values  of  [n*7r*  +  K  —  mrQ)  or  [  (n  +  l)*Tr*  +K  —  (n  +  l)irQ] 
that  is,  only  after  the  divergence  speed  has  been  exceeded.  Hence, 
as  in  the  infinite  panel  case,  the  finite  panel  in  inccxnpressible 
flow  yields  divergence  below  flutter.  The  above  result  has  also 
been  shown  more  generally  for  an  infinite  number  of  inodes  by 
Flax  in  Ref. 6. 

For  the  present  panel,  the  four  roots  T  associated 
with  coupling  between  the  8*"^  and  9^^  modes  were  computed  using 
Eqs.53  and  50  for  various  values  of  Q.  The  four  corresponding 
dimensioned  roots  s  are  shown  in  Fig. 23  for  various  speeds.  At 
low  speeds,  these  roots  represent  two  undamped  steady  sinusoidal 
oscillations.  Then  static  divergence  sets  in  for  a  very  narrow 
band  of  speeds  from  101.55  to  101.77  mph.  Above  this  band, 
two  undamped  steady  sinusoidal  oscillations  are  again  present. 
Finally,  at  speeds  above  103.13  mph,  flutter  instability  occurs. 

The  rate  of  growth  associated  with  these  instabili¬ 
ties  is  of  interest.  For  this  finite  panel,  the  amplification 
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ratio  A  is  Introduced  which  is  defined  now  as  the  ratio  of 
amplitudes  during  a  time  Interval  of  27r/21.3  -  0.295  sec. 
This  time  interval  represents  the  time  for  one  oscillation 
at  the  flutter  frequency  and  is  chosen  in  an  attempt  to 
make  the  A  of  the  finite  panel  somewhat  comparable  to  the  A 
of  the  infinite  panel.  This  amplification  A  is  then, 

A  =  e  (56) 

and  is  shown  plotted  versus  speed  U  in  Fig. 24.  Again  it  is 
seen  that  the  divergence  Instability  is  a  very  mild  one 
compared  to  the  strong  instability  associated  with  the 
flutter  behavior.  This  fact  together  with  the  narrow  speed 
band  associated  with  this  divergence  seems  to  indicate  that 
flutter  will  be  the  dominant  instability,  at  least  for  this 
simplified  two>mode  example. 

The  steady  oscillations  occurring  both  below  and 
above  the  divergence  speed  bands  were  investigated  in  more 
detail.  Figure  25  presents  a  plot  of  the  frequencies  co 
versus  air  speed  as  obtained  from  the  basic  roots  s  of  Fig. 2 
The  mode  shapes  associated  with  these  frequencies  can  be  ob¬ 
tained  from  Eqs.48.  Using  the  first  of  these  equations,  it 
can  be  shown  that  for  a  harmonic  oscillation  of  the  form, 

4)t  (57) 

The  associated  mode  q9(t)  is  given  by 

(58) 

where. 
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8n 


IT 

iQf 


K-nJTQ) 


(59) 


and  the  frequency  <0  is  the  absolute  magnitude  of  a  pure 
imaginary  pair  of  roots  s  such  as  occur  in  Fig. 23.  The 
nondimensional  root  s  appearing  above  is  related  to  s 
through  Eq. 47.  Then,  the  deflection  shape  for  these  steady 
oscillations  is  given  by. 


Thus,  the  two  standing  wave  mode  shapes  sin  mrx/i  and  sin 
(n  +  l)Trx/£  are  coupled  together  90  degrees  out  of  phase  by 
the  aerodynamic  damping  forces.  This  is  true  for  all  the 
steady  oscillations  occurring  below  or  above  the  divergence 
speed  band. 

The  mode  coupling  factor  fa  plays  an  interesting 
role  in  Eq.60.  Values  of  this  factor  fa  *  0  and  fa  ■  co 
represent  pure  standing  waves  in  the  sin  tarx/l  and  sin 
(n  +  l)7rx/£  modes,  respectively,  while  values  of  fa  ■  +  1 
and  fa  ■  —  1  represent  approximately  travelling  waves  in 
the  upstream  and  downstream  directions,  respectively.  This 
latter  fact  can  be  seen  by  rewriting  the  deflection  shape 
Eq.60  equivalently  as. 


Wfa.t)  =  04  lot  _  cea(h^i)'W 

+  cot  l| 

(61) 
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Hear  the  center  portion  of  a  long  panel,  . 

OircL  g  ^  7^71  (62) 

This  approximation  is  particularly  good  for  long  panels 
where  n  is  large.  Then,  Eq.61  reduces  to 


W  I^iot  —  45*J 

W(a,t)  ^  [wt  -  +45*] 


for  =  f  i 

(63) 

For  f,_  =  -  1 


These  are  readily  identified  as  waves  travelling  to  the  left 
(upstream)  and  to  the  right  (dovmstream) ,  respectively. 

A  plot  of  the  factor  fa  versus  air  speed  for  this 
panel  is  given  In  Fig. 26.  The  behavior  of  these  steady  os¬ 
cillations  then  becomes  clear.  At  very  low  speeds,  they 
represent  more  or  less  pure  standing  waves.  As  the  speed 
increases,  they  rapidly  develop  into  two  travelling  waves, 
the  lower  frequency  branch  L  moving  upstream,  the  higher 
one  H  moving  downsteam.  At  the  divergence  speed,  the 
lower  frequency  branch  L,  drops  sharply  to  zero  indicating 
a  divergence  in  the  lower  mode  8  again.  At  the  upper  diver¬ 
gence  speed,  the  divergence  occurs  in  the  higher  mode  9. 
Above  the  narrow  divergence  speed  gap,  the  lower  frequency 
branch  LL  reappears  and  rapidly  becomes  a  travelling  wave 
now  moving  downstream.  At  the  flutter  speed,  the  two  travel¬ 
ling  waves  coalesce  and  give  rise  to  unstable  oscillations. 

A  sketch  of  the  oscillations  at  flutter  is  given  in  Fig. 27. 
The  travelling  wave  character  of  these  oscillations  is 
clearly  apparent. 
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The  characteristics  at  flutter  are  calculated 
fron  this  tv/o  degrees  of  freedom  analysis  as. 

Flutter  speed  Up  *  103.1  mph 

Flutter  frequency  -  21.3  rad/sec  -  3.4  cps 

Wavelength  X  =  2i/(n  +  -j)  ■  24.5  in  (65) 

Wave  speed  c  4*  U  -  uiS/lir  »  6.9  ft/scc 
These  are  to  be  compared  with  the  corresponding  values  for 
the  infinite  panel  given  by  Eq.23.  The  flutter  speeds  are 
similar;  however,  the  flutter  frequency  and  wave  speed  are 
somewhat  lower  here. 

The  general  behavior  of  the  two  degrees  of  freedom 
finite  panel  analysis  closely  resembles  the  infinite  panel 
analysis  described  in  Section  3.2.  The  divergence  speed, 
the  mild  character  of  divergence,  the  undamped  travelling 
waves  moving  in  opposite  directions  below  the  divergence 
speed  and  in  the  same  direction  above  the  divergence  speed 
(compare  with  Fig. 17),  and  finally  their  merging  together  to 
form  a  strong  flutter,  arc  all  present  here.  The  character¬ 
istics  at  divergence  and  flutter  are  also  quite  similar. 
Presumably,  the  presence  of  more  modes  and  some  additional 
internal  damping  would  bring  the  results  of  the  finite  and 
infinite  panel  analyses  even  more  closely  in  line  with  one 
another,  particularly  for  instabilities  occurring  at  high 
values  of  n.  Thus,  a  smooth  transition  between  the  finite 
panel,  standing  wave  analysis  at  low  values  of  n  to  the  in¬ 
finite  panel,  travelling  wave  analysis  at  high  values  of  n  is 
observed. 

To  conclude  this  section,  it  might  be  appropriate 
to  give  here  a  simple  criterion  for  these  finite  panels  on 
clastic  foundations.  Using  Eq.51,  the  onset  of  divergence 
can  be  defined  by  the  nondimensional  parameter  Q  which  is  a 
function  only  of  K.  This  is  shown  plotted  in  Fig. 28.  For 
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values  of  K  greater  than  2000,  the  curve  shovm  there  can 
be  approximated  by 

q  «  I.7SS  k’'* 

n  ^  .Z+Z  ^  (nearest  (67) 

integer) 

The  criterion  Eq.66  reduces  identically  to  that  given  by 
Eq.AO'for  the  infinite  panel.  As  discussed  in  Sections  J.2 
and  3-3,  these  divergence  criterions  represent  only  a  lower 
bound  on  the  occurrence  of  Instability.  The  actual  in¬ 
stability  is  more  apt  to  be  Che  stronger  flutter- type  in¬ 
stability  which  occurs  somex^hat  above  the  divergence  speed 
quoted  here. 
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4.  CaiPARISOM  OF  ^PERIMENT  AIJD  THEORY 

The  experimental  test  results  obtained  for  this 
panel  were  described  in  detail  in  Section  2.5,  and  were 
summarized  in  Section  2.5f.  These  results  indicated  definite 
flutter  with  the  following  approximate  characteristics. 

Up  ^  95  mph 

\  ^  27  inches 

c  +  U  35  f c/sec 

a>  14.4  cps. 

The  waves  generally  appeared  to  be  travelling  downstream, 
although  there  seemed  to  be  some  standing  wave  component 
present  near  the  center  of  the  panel.  Below  the  occurrence 
of  flutter,  the  panel  assumed  a  small  static  sinusoidal  de¬ 
flection  shape  of  about  two-foot  wavelength  over  the  rear 
half  of  the  panel. 

The  theoretical  results  for  this  panel  were  de¬ 
scribed  in  detail  in  Section  3,  and  showed  essentially  the 
same  behavior  from  both  the  infinite  panel  and  finite  panel 
point  of  view.  Using  the  infinite  panel  analysis,  these  re¬ 
sults  indicated  a  strong  travelling  wave-type  flutter  de¬ 
veloping  with  the  following  characteristics. 

Up  »  105.5  mph 

^  «  25  inches 

c  +  U  “  11.4  ft/sec 

u)  ■  5.5  cps 

In  addition,  a  much  milder  Instability  of  static  divergence 
type  was  also  shown  to  exist  for  this  lightly  damped  panel 


at  l<3wcr  speed  of  101.5  mph  and  with  essentially 

the  same  X  "  25  inches. 

Comparing  these  experimental  and  theoretical  re¬ 
sults,  it  is  seen  that  there  is  good  agreement  in  the  flutter 
speed  and  wavelength,  but  rather  poor  agreement  in  the  wave 
speed  and  frequency  at  flutter.  Below  the  flutter  speed, 
the  observed  small  static  deflections  agree  well  in  wave¬ 
length  with  those  of  theory. 

The  poor  agreement  found  here  in  wave  speed  and 
frequency  at  flutter  is  typical  of  other  aeroelastic  phe¬ 
nomena  in  which  the  flutter  frequency  is  less  accurately 
predicted  than  the  flutter  speed.  Looking  at  the  variation 
of  wave  speed  c  +  U  with  air  speed  given  in  Fig.  17  or  at  the 
variation  in  frequency  to  with  air  speed  given  in  Fig. 25,  it 
is  seen  that  very  small  changes  in  air  speed  have  great  ef¬ 
fects  on  the  wave  speed  and  frequency  in  the  vicinity  of  the 
flutter  speed.  The  prediction  of  c  +  U  and  u>  seems  here 
generally  more  difficult  than  that  of  air  speed. 

Some  specific  factors  which  may  have  contributed 
to  the  poor  agreement  in  wave  speed  and  frequency  are  given 
below. 

a)  The  static  depression  at  the  leading  edge. 

This  shortened  the  effective  length  of  the  panel  and  may 
have  introduced  midplane  tension  and  structural  nonlinear i- 
tics  into  the  panel  behavior.  It  is  felt  that  this  effect 
was  the  primary  reason  for  the  scatter  in  the  experimental 
results  at  high,  medium,  and  low  angles  of  attack.  Also, 
this  probably  contributed  to  the  appearance  of  the  large 
static  deflections  of  Runs  114  through  121.  More  care 
should  be  taken  to  maintaining  uniform  pressure  and  velocity 
distributions  over  the  panel  surface  by  better  wall  con¬ 
touring  and  venting. 
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b)  The  cavity  effect  of  the  still  air  on  the 
bottom  side  of  the  panel.  The  effect  of -still  air  of  in¬ 
finite  depth  on  the  bottom  side  of  the  panel  was  investi¬ 
gated  in  Appendix  A  and  found  to  be  negligible  for  this 
panel.  Hov;ever,  it  nay  bo  that  for  the  finite  depth  of 
about  2  1/2  inches  here,  the  cavity  resonance  effect  may 
be  more  appreciable.  Also,  the  large  amplitudes  of  the 
flutter  deflections  may  have  had  an  influence  here. 

c)  Buffeting  at  the  trailing  edge.  There  nay 
have  been  some  forcing  oscillations  caused  by  rough  flow 
at  the  trailing  edge,  particularly  when  the  edge  was  free. 
The  somev/hat  similar  experimental  results  obtained  for  the 
case  of  free  an4-  pinned  edges  seen  to  minimize  this  source 
of  disturbance. 

d)  Boundary- layer  effects.  The  effects  of  a 
boundary  layer  in  influencing  panel  flutter  behavior 
have  recently  been  suggested.  Some  investigations  have 
been  carried  out  by  Miles  in  Ref. 7,  and  more  recently, 
Landahl  in  Ref. A,  studied  the  stability  of  a  boundary 
layer  in  the  presence  of  a  flexible  wall.  .It  is  believed 
this  specific  panel  lies  outside  the  range  of  any  such  sig¬ 
nificant  effects. 
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5.  CONCLUSIONS 

The  present  report  has  investigated  the  subsonic 
flutter  of  panels  on  continuous  elastic  foundations.  In¬ 
teresting  phenomena  have  been  observed  both  experimentally 
and  theoretically. 

Experimentally,  definite  flutter  was  observed  for 
a  specific  panel  in  this  subsonic  flow.  The  flutter  was  of 
a  travel ling-wave  type,  although  there  seemed  to  bo  some 
standing  wave  component  present  near  the  center  of  the 
panel.  Below  the  occurrence  of  flutter,  the  panel  assumed 
a  small  static  sinusoidal  deflection  shape  over  the  rear 
half  of  the  panel. 

Theoretically,  it  has  been  shown  that  a  finite 
panel  analysis  gives  essentially  the  same  behavior  as  the 
infinite  panel  analysis  for  high  values  of  n  (or  f/^) • 

Thus,  a  smooth  transition  is  demonstrated  between  the 
finite  panel,  standing  wave  analysis  at  low  values  of  n  to 
the  infinite  panel,  travelling  wave  analysis  at  high  values 
of  n.  It  has  also  been  demonstrated  that  although  a  mild 
static-divergence-type  instability  may  exist  for  these 
panels,  the  dominant  instability  will  be  a  travelling  wave, 
flutter-type  instability.  The  role  of  additional  internal 
damping  was  also  investigated  for  these  panels,  and  was 
shown  to  have  an  interesting  destabilizing  effect  in  the 
region  between  the  divergence  speed  and  the  flutter  speed. 

Comparison  of  the  experimental  results  with 
theory  showed  good  agreement  in  the  prediction  of  a  travel¬ 
ling  wave- type  flutter,  the  flutter  speed,  and  the  wave¬ 
length.  However,  there  was  rather  poor  agreement  in  the 
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prediction  of  the  wave  speed  and  frequency  at  flutter.  Out¬ 
side  of  the  general  difficulty  of  predicting  the  wave  speed 
and  frequency  as  accurately  as  the  flutter  speedy  this  dis¬ 
crepancy  was  attributed  possibly  to  limitations  in  the  ex¬ 
periment  arising  from  a  static  depression  at  the  leading 
edge  of  the  panel,  the  cavity  effect  of  the  still  air  on 
the  bottom  side  of  the  panel,  buffeting  at  the  trailing 
edge,  and  boundary- layer  effects. 

Much  further  work  remains  to  be  done  on  this 
general  problem  area  to  understand  it  more  fully.  More 
careful  experiments,  better  agreement  in  wave  speeds  and 
frequencies,  experimental  verification  of  the  interesting 
relationship  of  divergence,  flutter,  and  internal  damping, 
the  effects  of  still-air  cavities,  the  effects  of  boundary 
layers,  all  need  yet  to  be  done.  The  related  problem  of 
subsonic  flutter  of  panels  without  elastic  foundations  is 
also  of  considerable  interest.  Some  interesting  work  has 
been  done  by  Greenspon  and  Goldman  in  Ref. 8  and  again  the 
existence  of  experimental  subsonic  flutter  was  revealed. 
These  subsonic  investigations  should  be  of  interest  in  the 
problems  of  hydroelastic  panel  flutter,  axially-symmetric 
cylinder  flutter  at  low  speeds,  and  panel  flutter  of  in¬ 
flatable  structures  at  low  speeds. 

The  extension  of  this  work  to  the  supersonic 
flutter  of  panels  on  continuous  clastic  foundations  would 
also  be  of  considerable  interest.  Here  again.  Miles'  in¬ 
finite  panel  theory  of  Ref.l  can  be  used  as  a  reference, 
and  the  effects  of  a  finite  panel,  internal  damping,  etc. 
can  again  be  investigated  and  compared  with  experiment. 
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APPENDIX  A 


INFLUEfICE  OF  STILL  AIR  OH  BOTTOf-I  SIDE  OF  PANEL 


In  the  theory  developed  in  Section  3,  the  pres¬ 
sure  acting  on  the  bottom  side  of  the  panel  was  assigned 
constant  at  the  value  p  oo  •  It  is  of  interest  to  investi¬ 
gate  the  additional  pressure  caused  by  the  presence  of 
still  air  on  this  bottom  side. 

a)  Infinite  Panel 

The  pressure  on  the  bottom  surface  of  a  sinu¬ 
soidal  wave  travelling  with  wave  speed  "c  through  still  air 
has  been  shown  by  Miles  in  Ref.l  to  be 


^  c  w.  e 


(A-l) 


In  the  above,  x  represents  the  distance  relative  to  the 
panel  and  “c  represents  the  wave  speed  relative  to  the 
panel.  To  convert  the  above  results  to  the  axis  system 
used  by  Miles  in  his  analysis  and  repeated  again  in 
Section  3.2,  one  notes  that 

X  “  X  —  Ut 

•  c  “  ^  —  U 


(A-2) 

(A-3) 


where  x  and  c  refer  to  a  fixed  axis  system,  and  the  plate 
itself  moves  with  velocity  U  in  the  negative  x  direction. 
Placing  Eqs.A-2  and  A-3  into  Eq.A-1  gives  the  additional 
pressure  from  the  bottom  side  of  the  panel  as 


(ot-x; 


(A- A' 
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The  above  pressure  is  added  to  that  for  the  top  side  of  the 
panel  given  by  2q.ll,  and  then  substituted  into  the  moving 
panel  equation,  Eq.9.  Upon  nondimensionalization  as  be¬ 
fore,  one  now  obtains  the  equation 

in  place  of  the  previous  Eq.l2. 

Solving  Eq.A-5  explicitly  for  one  obtains  the 
following  expressions  for  the  wave  speed  relative  to  still 
air,  c, 

±  ]  (A-6) 


and  for  the  wave  speed  relative  to  the  panel  c  +  U  •  c. 


These  expressions  arc  to  be  compared  with  the  previous  ex¬ 
pressions  Eqs.l7  and  18,  which  neglected  the  air  pressure 
variation  on  the  bottom  side  of  the  panel. 

It  is  readily  seen  from  the  above,  that  negative 
imaginary  values  of  c,  hence  instability,  occur  now  for 


H  ^  I  » ^  ' 

c,  4 


(A-8) 


This  criterion  is  indicated  in  Fig.lA  by  the  dashed  lines. 
For  values  of  p.  <  0.3,  this  new  flutter  criterion  is  the 
same  as  the  previous  one  given  by  Eq.l9.  At  the  higher 
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values  of  p..  It  approaches  1.414  times  the  divergence  cri¬ 
terion.  The  corresponding  wave  speed  relative  to  the  panel 
at  flutter  is 


O  -K  IT 


(A-9) 


For  the  present  panel,  the  characteristics  at 
flutter  with  still  air  on  the  bottom  side  are  calculated 
as 


Up  -  105.3  mph 
X  -  25  in 


c  +  U  -  10.5  ft/sec  (A- 10) 

to  “  ^  (c  +  U)  -  31.7  rad/sec  **5.1  cps 
IL  -  0.079 

Comparing  these  with  the  corresponding  ones  £q.23,  it  is 
seen  that  there  is  a  negligible  change  in  flutter  speed, 
and  only  a  slight  modification  to  the  wave  speed  and  fre¬ 
quency. 

The  divergence  characteristics  remain  unaltered 
by  the  presence  of  this  still  air.  Also,  the  general  char¬ 
acter  of  transition  between  the  divergence  speed  and  flutter 
speed  is  identical  for  the  low  p.  of  this  panel.  At  the 
higher  p's,  some  modification  might  be  present. 

The  theory  of  this  appendix  would  be  important 
for  panels  with  high  p's  such  as  might  result  for  very 
light  panels,  for  very  soft  spring  foundations  K  (hence 
large  heavy  fluid  media  such  as  water,  where 

the  fluid  was  present  both  on  the  top  and  bottom  side  of 
the  panel.  For  cases  where  water  is  present  only  on  the  top 
side  and  air  is  on  the  bottom  side,  the  theory  of  Section 
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3.2  would  be  more  applicable. 

The  analysis  presented  here  considered  the  still 
fluid  to  bo  large  in  depth.  The  effect  of  a  small  depth 
cavity  on  the  pressure  distribution  was  not  investigated 
hero.  Conceivably,  it  could  give  a  greater  influence  on 
the  flutter  characteristics. 

-•b)  Finite  Panel 

The  additional  pressure  from  the  vibration  of  the 
bottom  side  of  the  panel  in  still  air  is 


(A-11) 


This  represents  purely  an  additional  virtual  mass  term. 
Carrying  through  the  analysis  as  in  Section  3.3  will  re¬ 
sult  in  the  set  of  nondimcnsional  differential  equations 


(A- 12) 

• 

The  only  difference  between  the  above  equation  and  the  previ¬ 
ous  Kq.46  is  the  presence  of  a  factor  2  in  the  d*q^dt  * 
term,  representing  a  doubling  of  the  virtual  mass.  The 
characteristic  equation  for  the  coupling  of  the  nth  and 
n  +  1th  nodes  is  given  again  by  Eq.49,  but  now  the  coeffici¬ 
ents  A«  and  Aa  would  reflect  this  doubled  virtual  mass. 

For  the  present  panel,  the  four  roots  ’s  were  again 
computed.  The  divergence  speed  gap  is  unaffected  and  remains 
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between  101.55  and  101.77  nph  as  before.  The  characteris¬ 
tics  at 'flutter  arc  now  computed  as 

Flutter  speed  Up  *»  103.07  mph 

Flutter  frequency  tOp  ■  19.3  rad/soc  -  3.2  cps 

Wavelength  X  •*  24.5  in 

Wave  speed  c  +  U  “  (dK/Itt  6.5  *ft/sec 

Comparing  these  with  the  corresponding  Eq.65,  it  is  again 
seen  that  there  is  a  negligible  change  in  flutter  speed 
and  only  a  slight  modification  to  the  flutter  frequency. 


-13) 
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Fig.  6  Flutter  Test  Log 
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Pressure  Relative  to  Atmospheric 
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Fig.  10  Summery  of  Flutter  Observations 
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Wp  cs  12.2  c.p.s. 
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Fig.  11  Ftlin  Log 
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Fig, 14  Flutter  and  Divergence  Criteria  (Infinite  Panel 
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Fig. 16 
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F1^.19  Amplification  for  Undamped  Panel  (Infinite  Panel) 


70 


AmpllficaClon 


< 

I 


Air  Speed  ~V  (M.P.H.) 


72 


Nonditnenalonal  Parameter 


Fig. 22  Diver 
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Fig.  23  Roota  vs.  Speed  (Finite  Panel) 
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